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Abstract 

The effect of excluded volume interactions on the structure of a polymer 
in shear flow is investigated by Brownian Dynamics simulations for chains 
with size 30 < < 300. The main results concern the structure factor 
5(q) of chains of N=300 Kuhn segments, observed at a reduced shear rate 
/3 = ■jr = 3.2, where 7 is the bare shear rate and r is the longest relaxation 
time of the chain. At low q, where anisotropic global deformation is probed, 
the chain form factor is shown to match the form factor of the continuous 
Rouse model under shear at the same reduced shear rate, computed here for 
the first time in a wide range of wave vectors. At high q, the chain structure 
factor evolves towards the isotropic equilibrium power law q~^^'^ typical of 
self-avoiding walk statistics. The matching between excluded volume and 
ideal chains at small q, and the excluded volume power law behavior at large 
q are observed for q orthogonal to the main elongation axis but not yet for 
q along the elongation direction itself, as a result of interferences with finite 
extensibility effects. Our simulations support the existence of anisotropic 
shear blobs for polymers in good solvent under shear flow for /3 > 1 provided 
chains are sufficiently long. 
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I. INTRODUCTION 



Excluded volume (EV) and hydrodynamics interactions (HI) have been recognized as 
fundamental physical phenomena to explain the peculiarities of single chain behavior in di- 
lute polymer solutions. For a solution at rest, their effects on the single chain structure 
and dynamics have been carefully investigated over the years both theoretically and exper- 
imentally fTj-Q. When a dilute solution is subjected to a flow field, the situation is much 
more complex. In addition to finite extensibility (FE) effects to be considered at high strain, 
the role of EV and HI effects is still largely debated. This situation certainly arises from 
the difficulties to solve the kinetic equations for the dynamics of the simple coarse-grained 
models when EV and HI are considered. Some progresses have been made concerning the 
HI effects in absence of EV interactions by the use of a self consistent preaveraging approx- 
imations 0, a non equilibrium version of the more famous equilibrium Zimm theory [Q. 
Renormalization Group Theory techniques have also been applied to solve suitable dynamic 
equations for the chain The main focus of those investigations was the rheological 

response of the system although some results on the chain structure where also provided 
|lO,P]. An attempt to introduce EV effects in the Rouse model has also been performed by 



perturbation theory expansion [|13 



On the experimental side the situation is similar in that measurements of single chain 
properties in flow are considerably more difficult than at rest. However, several scattering 
experiments have been performed, mainly Light Scattering (LS), with the aim to investigate 
the global properties of the chain such as the orientation and the deformation of the coil by 
the flow ||l4|-pO[] . The effect of EV interaction over global chain properties has been discussed 



in a series of recent LS studies |18-20| where some systematic trends were reported. For 



dilute solutions in shear flow at a flxed reduced shear rate, a systematic weakening of the 
alignment of the main axis of the chain in the flow direction has been observed for increasing 
solvent quality. At a given reduced shear rate, the measured expansion ratio indicates that 
good solvent chains are less deformed than theta chains Inspection of the flow effects on 
the internal structure of the coil can only be obtained by Small Angle Neutron Scattering 
(SANS) and indeed some data have been already published |]T3[. However experimental 



limitations are here even more severe than for LS, and a quantitative study for various 
solvent qualities is still missing. 

New experimental insight in the single chain behavior under flows and stretching forces 
has been provided recently by fluorescence techniques on biological molecules like DNA 
pl|-ROl. Two papers on the single chain dynamics under shear flow have been published 



recently by S. Chu and coworkers p9[ and by D. Wirtz and coworkers in both cases 
employing DNA chains. In the former paper,the mean and time fluctuations of the fractional 
extension of DNA under shear flow in the flow direction was reported for reduced shear rates 
larger than one and culminating at /3 = 76. In the second paper, the authors have shown 
that the dynamics of a DNA chain under weak shear flow {f3 ~ 0.1) presents an unexpected 
wide range of relaxation times related to the particular initial conflguration prior to the flow 
inset. 

An interesting alternative to the experimental tools in the context of polymer physics 
are the numerical techniques such as Monte Carlo (MC), Molecular Dynamics (MD) and 
Brownian Dynamics (BD) simulations. They can be used to solve simple models for the 
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system of interest and therefore test or suggest theoretical developments or stimulate new 
experiments. On the other hand simulations can help interpreting experimental data in a 
more complete and coherent fashion. 

Many numerical studies of the behavior of a single chain under shear flow have already 
appeared by various techniques such as BD pT|^6|, MD [p7|-^ and even non-equilibrium 



MC [^^ . To mention a recent simulation work |3^, BD studies for a chain without HI 
and EV found a reasonable agreement with the infinite chain Rouse model predictions but 
came to the disturbing conclusion that, while the relative elongation of DNA chain in the 
flow direction is correctly explained by the chain model, the same model fails to explain 
the global expansion ratio of a polystyrene chain which is experimentally much lower than 
predicted. 

To our knowledge, the role of the excluded volume effect on the structure of the chain 
under SF has not been systematically studied by simulations. In the present paper we try to 
elucidate this effect by comparing the analytical results of the Rouse model (in the continuum 
limit) with BD data for a suitable microscopic model of flexible chain with excluded volume 
interactions and unavoidable FE features. 

The structure of a single chain in shear flow is often discussed by exploiting an analogy 
with the well known structure of a single chain subjected to a pair of equal and opposite 
stretching forces applied to its ends. In this case the " macroscopic" behavior of long chains, 
as measured by the elasticity law, is characterized by a linear hookean regime for weak forces 
followed by a non linear power law regime for strong forces and ultimately by specific finite 
extensibility effects. The crossover between the first two regimes is linked to the appearance 



of the so called "tensile" blobs in the chain structure P3|,3. 

A very clear signature of the blobs can be detected in the chain structure factor as a 
crossover from ideal behavior at small q's to excluded volume statistics at large q's. This 
has been shown long time ago by Neutron Scattering experiments for "thermal" blobs in the 
semi-diluted solutions at rest ^M, and recently by MC simulation for "tensile" blobs in the 



stretched chain problem In the latter case, the crossover is detected for q orthogonal 
to the external force and can be ascribed to the fact that scattering at small q's in these 
directions comes from pairs of monomers which are far apart in the elongation direction and 
therefore do not interact directly through the EV potential. 

For the chain in shear flow Onuki has formulated a blob model |^6[ as a simple extension 



of the Pincus-deGennes theory for the stretched chain case P3|,3. The blob size is related 



to the bare shear rate and the chain at high shear rates is imagined as a string of shear 
blobs aligned most of the time in the flow direction. Such a structure which corresponds to 
an aligned fluctuating rod, should yield a scattering signal for q in the flow direction (flow 
and elongation directions provide essentially the same scattering signal at high enough shear 



rate |3^) similar to the stretched chain signal for q aligned in the stretching direction ||45 



For the stretched chain such behavior is described quantitatively by the Rouse model with 
a suitable choice of the chain elongation and the longitudinal and transverse fluctuations 
,57]. In the present paper we test these ideas on the form factor of a single EV chain 



in flow for shear rates which are sufficiently high to reach the plausible shear blob regime 
but also sufficiently small to remain in the scaling regime, i.e. avoiding finite chain effects. 
Our strategy is to analyze to which extend the low q behavior resembles the single chain 
structure factor of ideal chains in the presence of shear flow, which can be derived exactly, 
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and whether the high q regime reproduces the structure of unperturbed EV chains. In 
order to observe such phenomena taking place on different length scales, we need sufficiently 
long chains and therefore we have performed BD simulations of EV chains in absence and 
presence of shear flow. 

The paper is organized as follows. In section |I| we set up the essential phenomenology 
concerning a chain in shear flow and define the relevant structural characteristics of the 
chain. Section describes the procedure to solve the Rouse model in flow and provides 
the results for the specific SF case. This is not a totally new analytical development as 
we are aware of at least two previous works on it [^,^. The novelty here is however 
the calculation of the structure factor in the whole range of q vectors and the use of such 
behavior to pinpoint the signature of shear blobs in the EV chain structure factor. In the 
next section |V| we describe the model for EV chains and the BD technique used. In section 
we discuss the results of BD simulations. Finally section is devoted to our conclusion 
and perspectives. 



II. PHENOMENOLOGICAL FRAMEWORK 

When a dilute polymer solution is subjected to a simple shear flow, characterized by the 
velocity field u = 7?/x where7 is the shear rate, the chains in the solution are oriented and 
deformed according to the flow. Those phenomena can be measured through the anisotropy 
arising in any tensorial quantity associated to the chain, such as for instance the gyration 
tensor G and the order parameter tensor O defined as 

G = ^ E < (R. - R,)(R^ - R,) > (1) 



i=l 



U 



where Rj are the coordinates of the i-th monomer in the chain, Uj = Rj+i — and 1 is the 
unit tensor. At equilibrium the system is isotropic and any tensorial quantity, say A reduces 
to a scalar quantity. Under shear flow, symmetry imposes A^z = = = A^y = 0. 
Therefore there are at most four independent components of any tensor to be monitored, 
namely the three diagonal elements and the off-diagonal A^y = Ay^. The orientational angle 
Xa, defined through the relation 

cot{2xA) = (3) 

measure the rotation around z of the principal axes (I, II, III) of the tensor A with respect 
to the flow axes (x,y,z). In shear flow, Axy starts linearly with 7 while the first contribution 
to Axx — Ayy is of order 7^. Therefore the linear (Newtonian) regime is characterized by 
Xa = Tc/4. Outside the linear regime xa decreases to zero for increasing shear rate. 
The deformation ratios, defined by 
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measures the amount of deformation of the chain in flow, either in the flow reference frame 
{a = X, y, z), or in the molecular reference frame {a = I, II, III). When defined in the flow 
reference frame, the first contribution to those quantities is of order 7^, so that the Newtonian 
regime is characterized by the absence of deformation. Outside the linear regime, the chain is 
elongated in the flow direction {6Axx > 0) and compressed in the in-plane gradient direction 
{SAyy < 0), while the neutral (out of plane) direction is only shghtly decreased {6Azz ~ 0). 



III. STRUCTURE FACTOR OF THE ROUSE CHAIN UNDER SHEAR FLOW 

To compute the structure factor of an ideal gaussian chain (Rouse chain) we adopt 
the continuous chain model introduced by Edwards in which the only energy term is a 
quadratic potential acting between nearest neighbors along the chain. For such a model the 
equilibrium probability distribution in configurational space is 

where x = SkBT/d"^, represents the positions of the chain monomers, < n < is 
the continuous monomer index and d is the single bond variance of the model. The steady 
state distribution in presence of a generic homogeneous solvent flow field can be obtained 
in terms of the normal modes of the model (Xp, (p = 0, 00)) 0,0. For homogeneous flows 
to which we limit our discussion, the normal modes are defined by the same transformation 
as at equilibrium. In those new coordinates the distribution is factorized over the normal 
modes 

00 

^({x},t) = n^p(Xp,t) (6) 

p=0 

In SF, the steady state distribution of the p-th mode is 



where the matrix jSp is defined as 



ft = l + ^(k + k-) + 2(^j'k.k'- (8) 

and k = Vu is the velocity gradient tensor. The derivation of eqs. (0) and (|) is given in 
appendix ^ 

The chain structure factor is defined as 

1 i-N i-N 

S(q, N,i) = — / dndm < e^'i-^— >^ (9) 
N Jo Jo 

where Unm = Rn — R-m and < ... >^ indicates a statistical averages over the non equilibrium 
distribution eq.(BD. In terms of normal modes we have 



5 



2E 

p=i 



pun pnm\ 
cos—- cos 



N 

and using the model solution eq.(^ we obtain, after a gaussian integration, 



(10) 



exp 



2kBT 



Xi 



-q- 



p=l 



cos- 



prm 
1^ 



cos- 



pnm 
1^ 



■ q 



(11) 



Substituting the expression of /3p from eq. (||) we get 



cxp ■ 



2kBT 
Xi 



S'o(n, m)g^ + -^Si{n, m)q • (k 



k^).q 



jS2{n,m)q - (k ■ k^) -q 
(12) 



with 



oo -j^ 

p=i t' 



cos- 



pim 

w 



cos- 



pirm 



1 2 



N 



A; = 0,1,2 



(13) 



In eq.(|T2|) the first term in the exponent on the r.h.s. is the equilibrium contribution leading 
to the well known Debye function. The other two terms arise from the presence of the flow. 
The three series above can be summed analytically as explained in appendix |B| and the 
result is 



So{a,b) 



TX 



\a — h\ 



5i(a,6) = ^(a-6)2 



:a + 6) 1 



\a — h\ 



(14) 
(15) 



^2(a,6) = ^(a-6)2{|a-6|^ + 5(a + 6) [(a + 6)(2 + + 6^) - Sa^ - 2a6 - 36^1 } (16) 

where we have defined reduced indices a = n/N^h = m/N and < (a, 6) < 1. 

Let us now introduce the reduced shear rate /3 as the product of the bare shear rate 
and the global relaxation time of the chain at equilibrium. To compare the theory with the 
experiments we define {3 as 



NAksT 



(17) 



where M is the molecular weight of the polymer, rjs the solvent shear viscosity, [r]] the 
intrinsic viscosity of the solution and the Avogadro number. In the Rouse model the 
intrinsic viscosity is [rf\ = Na/Mtis {N'^d?^l?>Q) so that 



7 N'dX 
ksT 36 



(18) 



which gives ri7 = 12(3/71 0. The chain structure factor in terms of the reduced shear rate 
is finally obtained as 
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N JO 
where f (a, b, j3) is the following tensor 

12(3 



da 



db e~ 



i{a,b,(3) = |a-6|l + 



lqq^:f(a,6,/3) 



12/5 



5i(a, 6) (A + A' ) + ^52(a, 6) (A ■ A^) 



(19) 



(20) 



and A = 7~^k. Expressions ([T9|) with (pOD are our main interest here. The double integral 
over monomers cannot be performed analytically and we resort to numerical method. Before 
presenting specific results we want to make some further considerations. Firstly we note that 
the tensor f entering in the structure factor depends on and 7 through (3 only. This implies 
the following universality for the structure factor 



N 



F{Rgq,(3) 



(21) 



where F is a universal scalar function depending on the geometry of the applied flow and 
Rg = JNd'^/Q is the equilibrium radius of gyration for the gaussian chain. Secondly we 



want to obtain the well known expression for the gyration tensor of the ideal chain under 



shear flow p8| , p2| from which we can derive the orientation and deformation of the chain to 
be used below. The small q expansion to order of the structure factor is directly related 
to the gyration tensor by 



When applied to eq.([T 

Nd"^ '•I 



N 

it provides 



1 - qq' : G + 0(g' 



G{N,f3) = / da / dbi{a,b,P) 

6 JO Jo 



Nd'^ 
~6~ 



-1 + —(A + A^ 
3 15^ 



H —(A ■ A^ 

105 ^ ' 



(22) 



(23) 



Note that at fixed (3 the gyration tensor has the same N dependence as at equilibrium (the 
scaling exponent is u = 0.5). The orientation of the chain, as given by eq.(^, when defined 
through the gyration tensor is 



cotg{2xg 



yy 



60/3 (3 



2G 



(24) 



xy 



105 1.75 
1.75 and /3-independent. This is the well 



implying that the orientational resistance nig 
known classical result. 

Other used quantities in characterizing chains under SF are the deformation ratios as 
defined in section 0. From eq.(||) we get 5G^^ = 16/35 = 0A57/3^;5Gyy = 5G^, = 
as is well known. Note that the above values of the deformation ratios in the laboratory 
reference frame imply that in the molecular reference frame both elongation and compression 
directions present deformations as follows 



^^'=35 



SGii = — 



35 



4/?2 + 7/3 




(25) 



(26) 
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with the following limiting behaviours 



SGi —(3^ 13 »- (27) 

35 ^4 ^ ^ 

49 7 
5Gn ^-^ = -0.35 (3»- (28) 

14 7 
= ^ -/3 /? « - (29) 

Now we go back to our original task, that is the calculation of the chain structure 
factor under shear flow. As already stated, the double integral over monomers in eq.(|19D 
with eg . (1201) cannot be performed analytically because of the exponential character of the 
integrand (to be compared with the double integral of the exponent which can be solved 
and provides the gyration tensor as in eq. (p3[) ). We applied standard Romberg's method 
[|50| to numerically evaluate those integrals in a wide range of g's and for given (3. In figure 
m we report the chain structure factor in its universal form {S{q){qRgY /3N vs qRg/^/S), for 
various values of (3 and for q along the elongation and the compression directions in the flow 
plane (molecular reference frame). We note as in the elongation direction the crossover to 
the high g's equilibrium power law behavior strongly depends on the shear rate. The Kratky 
plot in the compression direction presents a /3 dependent overshoot due to the compression of 
the corresponding gyration tensor component. This effects however saturates very soon (the 
curve for f3 = 3.2 is already very close to the behavior for (3 = 10 which is indistinguishable 
from the one at /? = 100) and the crossover to the high g's power law behavior is in practice 
P independent. 



In comparison with the stretched ideal chain |^ we note the absence of oscillations in the 
elongation direction of the sheared chain signal. Such oscillations are characteristics of an 
aligned rod with uniform density of scattering centers. Our present results show definitively 
that, even at very high shear rates, the sheared chain cannot be thought in terms of such 



simplified model at variance with Onuki scheme [^] . 

Another important characteristics of the structure factor of the stretched chain, related 
to the presence of tensile blobs, is the Lorentzian shape of the scattering function at in- 
termediate g's as discussed by Benoit et al. [47]. In this case one can easily show that, 
for large g, the scattering function for q along the stretching direction follows the equilib- 
rium Ornstein-Zernike behavior S{q) = 1/(1 -|- g^i?^/2) if plotted in terms of the effective 
wave vector q"^ = q"^ + 4/^^, where C,t is the tensile blob size. This led Pincus to for- 
mulate his scaling for stretched EV chains just changing the scaling exponent from the 
ideal to the EV value |43,45|. Under shear flow the chain structure factor at intermediate 



g's in the elongation direction has not a Lorentzian shape, but rather presents a behavior 
S{x)/N = 1/(1 + + a{(3)x'^^'^'>) for x = qRg/^{?,) > 1, in analogy with the scattering 
signal from near critical fluids under shear The fltting function a{f3) goes linearly from 
to 73 in the interval < f3 < 100, while a{f3) goes almost linearly in log{j3) from 0.62 
to 0.88 in the interval 3.2 < (3 < 100. Whether this form for S{qi) has some theoretical 
foundation remains unclear to us. 
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IV. EXCLUDED VOLUME CHAIN MODEL AND BROWNIAN DYNAMICS 

ALGORITHM 



The ideal chain model of previous section is only the starting point for a more realistic 
description of polymeric chains. When the aim is the study of dilute solutions at least 
two additional effects need to be taken into account namely the excluded volume and the 
hydrodynamic interactions between monomers. In this section we describe the model we 
have used to represent the excluded volume effect. In our previous work on stretched chains 



45| we have used a model of linear chains with rigid bonds and no angle hindrance. The 
excluded volume interaction where modeled by hard sphere interactions. The diameter a of 
each monomer was chosen as cr = 0.65 in units of the bond length. This model exhibits good 



solvent scaling |^5[ and is very suitable for Monte Carlo calculations. Here we deal with 
non equilibrium systems for which the stationary distribution under shear is unknown and 
therefore we cannot apply the Monte Carlo sampling procedure. The only way to perform the 
calculation is by a dynamical technique such as Brownian Dynamics(BD). For this technique, 
the model of the previous work is not very appropriate because of the strong discontinuity in 
the interactions. We thus replaced the hard sphere by a Lennard- Jones potential truncated 
at the minimum and shifted in such a way to have a continuous potential 

VEv{r) = VLj{r) - VlAtc) r < = (30) 

= r > (31) 

where vlj{t) = 4e[(o"/r)^^ — (a/r)^]. In our calculation e was the energy unit and a was 
chosen to be 0.65 as in the previous work. Moreover the imposition of the bond constraints 
in BD requires an iterative technique such as SHAKE [^] that for long chains is quite 
demanding in terms of computer time. We therefore prefer to use flexible rather then rigid 
bonds and represent the bond interaction by an harmonic potential with force constant x 
and minimum position d: Vhondir) = X (kl ~ d)"^ /2. Lennard- Jones interactions between 
nearest neighbors along the chain were not considered. 
The equation of motion of the i-th monomer is 

eR, = F, + f, + ek-R, (32) 

where Fj is the total force on monomer i deriving from the potential energy, and the prop- 
erties of the noise fj are 

< f^(t) > =0 (33) 
< f„(t)f,(t') > = 2kBT5{t - t')5,,l (34) 

The above equations can be integrated numerically by the simple finite differences scheme 
R,(t + h) = [l + Kk] ■ R,(t) + jYiit) + w,(t) + 0{h"^) (35) 

where the moments of the white noise Wj are 

< Wi(t) > = (36) 
< Wi(t)wj(t') > = 2hDQ5{t - t')5ijl (37) 
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and Dq = ksT is the single monomer diffusion coefficient. Higher order schemes has been 
proposed in the hterature which allow to take a larger time step but require more operations 
per time step. A further advantage of the simple first order scheme over higher order ones, 
is the possibility to use a uniformly rather than normally distributed noise, provided that 
the zero and second moment are correctly chosen. Indeed, it is easy to show that the error 
introduced by such a substitution is of order h^/'^ or higher, that is beyond the precision 
of the scheme itself. Use of uniformly distributed noise rather than a gaussian noise saves 
about 50% of the computer time. In absence of a systematic comparison on the overall 
efficiency of integration schemes we have preferred to use the simple one above. 

In all our calculations we have fixed fc^T = 1 and ^ = 1 which lead to Z^o = 1- The bond 
force constant has been chosen to be x = 40 and in order to obtain a unitary average bond 
length we had to fix d = 0.945. 

The presence of stiff springs, ensuring the rigidity of the bonds, and strongly repulsive 
core interactions, ensuring the self avoidance and avoiding the self crossing of the chain, 
impose a quite small time step for the stability of the numerical scheme. On the other hand 
the long relaxation time of poljTiiers and its rapid increase with the number of monomers 
{tn ~ A^^'^ for excluded volume chains without HI) would require a time step as large as 
possible. This is the main limitation of any dynamical scheme when dealing with polymers. 
In the present work we have chosen a time step h = 0.00025 in reduced units. We have tested 
in some cases that halving the time step does not change the results for the chain structure. 
Even with such small time step, numerical instabilities occasionally show up. These occur 
when two particles end a time step at a reciprocal distance considerably smaller then a. 
The following time step is then driven by very large repulsive forces which stretch the bonds 
related to the pair of particle at a distance considerably larger then the minimum of the 
bond potential and so forth. To cure this occasional pathology we checked at the end of 
any time step the minimum distance between any pair of particles and if it is found to be 
smaller then 0.43cr we reject the time step and we proceed the time integration for a time 
interval h in many smaller time steps. When the time interval h is reached the original time 
step is restored. With the chosen value for the original time step we have found that such 
pathological events happen on average any 2x10^ time steps and that dividing the time step 
by ten is always enough to escape from the instability. This is a craft made version of more 
sophisticated adaptive time step integration schemes which are however more demanding in 
terms of computer time. 

To further minimize the computer time per step we have used a Verlet neighbors list |5B 
which is very effective with such short range interactions. 



V. SIMULATION RESULTS 

A. Equilibrium scaling of the model 

Before proceeding to study the chain model under shear flow we need to characterize 
its equilibrium behavior. In particular we need to know at which length the chains start 
following the known static and dynamic scalings and we need to determine the prefactor in 
the scaling laws in order to fix the values of the bare shear rate for a given chain length 
to work at fixed f3. In figure |^ we report the radius of gyration Rg and the end-to-end 
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distance R versus the number of links in the chain (A^ — 1). The observed scahngs are 
Rg = 0.40448(A^ - 1)° ^ and R = im{N - 1)°-^ if we exclude N = 9,20. Equilibrium 
structure factors for the various chain lengths are shown in figure |^. In figure we report 
the universal plot {qRgYS{q)/N versus (qRg)- We also add the Debye curve for the ideal 
chain. We observe that, for > 30, the EV chains follow indeed a universal q^^^ power law 
as predicted on the basis of the static scaling and dimensional arguments 0. 

By construction, the chain centre of mass diffusion coefficient is Dcm = Dq/N oc N'^ 
where v' is the dynamic exponent, which provides a characteristic time of the chain tat = 
Rl/{6Dc^) = (0.40448)2(A^ - l^-^N/QDo ~ N'^^+\ Here u = 3/5 is the Flory scaling 
exponent for EV chains, and z/' = 1 as HI are absent. In order to check the consistency 
between global and internal dynamics we extracted the characteristic time of the latter by 
inverting the chain dynamic structure factor J(g, t) = S{q, t) / S{q) as we have previously done 
for MD chains |]5^ , |39| . Indeed scaling and dimensional arguments predict that I{t,q,N) be 
a universal function of {tq~^) with x = 2 + u'/u = 11/3 in the present case. Figure ^ shows 
the behavior of (tq^^^^) vs q for two values of / and for various chain lengths. We observe 
that for = 20 the chain does not exhibit the correct scaling of internal times while from 
= 30 on, both static and dynamics are in the scaling laws regime. The data we show are 
results of quite long runs as reported in table |. 



B. Chain structure at fixed /? 

In this subsection we describe the results for the chain structure under shear flow for 
various chain lengths at the same value of the reduced shear rate j3 that we defined as 
P = 'yRg/{6Dcm)- We have hmited our study to two values of /5, namely /3 = 3.2 as in our 
previous MD study and in several experiments [|1^,^], and f3 = 10. In tables ^ and \[T\\ , we 



report various technical details of our simulations and the results for the orientation and the 
deformation. We studied chains up to A^ = 300 which require very long runs (up to 3 10^ 
time steps for the EV case, see table fT|). 

As in our previous MD study of short chains (up to A^ = 50), we analyze the results 
in the molecular reference frame on account of the observation that the orientational angle 
depends on A^ and 7 through P only. This is also suggested by the Rouse model results 
of section |ITI| and by its extension when HI at the level of equilibrium preaveraging are 
considered Our new results for EV chains at /5 = 3.2 confirm this expectation within 
error bars, although short chains seems to be slightly less oriented as shown in figure |^. 
On the same figure we report the extinction angle related to the birefringence which also is 
found to be the same within error bars for all A^ at fixed f3. 

Before discussing the results for our chain model with EV interactions it is instructive to 
consider the same coarse grained model with the EV interactions switched off (e = in the 
LJ potential). Comparison of the chain structure factor of such a model with the continuous 
ideal model result eq.(0) (not shown here) validates the theory itself and the numerical 
procedure applied to compute the double integral. Moreover data for the orientation and 
the deformation of this model can be useful in understanding qualitatively the role of the 
finite extensibility effects due to the stiff harmonic potential which keeps the bond length 
nearly constant. Data for the orientation and deformations of ideal chains with A^ = 100 
and A^ = 300 at /3 = 3.2 and (3 = 10 are reported in table |T| and compared to the Rouse 
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model predictions. At (3 = 3.2, the orientation for all chain lengths is in agreement with 
the long chain limit while the deformations for N = 100 deviate considerably from the ones 
for N = 300 (except in the direction II) which are instead in agreement with the Rouse 
prediction. Interestingly, the main eigenvalue Gj of the gyration tensor follows an apparent 
power law A^^'^ qualitatively similar to the behavior observed for MD data on shorter chains 
with EV and HI At /5 = 10 both chains exhibit finite extensibility effects even on the 
orientation. 

Now we turn to the EV results. In table |Tl| we collect the data for the orientation and 
the deformation of our chains. We also report data for the longest chains studied in ref. 



by molecular dynamics. We note at /5 = 3.2, that the 300-bead chain appears to be slightly 
more oriented and deformed than the shorter chains, although the noise is quite large. At 
(3 = 10 instead, the 300-bead chain is considerably more oriented and deformed than the 
100-bead chain, showing an important influence of the finite extensibility. Comparison with 
previous MD data for N=50 shows a general agreement for all quantities. 

From these results and those of table ||, it appears reasonable to assume that the N=300 
results at /3 = 3.2 are representative of the global properties of long chains, namely that 
the remaining small systematic finite chain effects are masked by statistical errors of a few 
percents (due to finite statistics). Comparing EV chains to 9 chains at the same reduced 
shear rate, we observe for EV chains a less marked alignment with respect to the flow (larger 
X value) and a smaller global deformation, in qualitative agreement with recent LS results 



An important question concerns the universal character of the structure factor under 
shear flow. It is well known that at equilibrium the quantity S{q,N)/N is a universal 
function of qRg only, for any coarse-grained model and for long enough chains 0. In the case 
of the chain stretched at both ends, we have shown that different universalities in terms of 
qRg hold separately in the longitudinal and the transverse directions |^. Alternatively one 
could investigate self-similarities in terms of q\/Ga where are the diagonal components 
of the gyration tensor. For the chain stretched by its ends, it is well known that fluctuations 
of the end-to-end distance (and therefore Ga) at fixed reduced force scale with as the 
equilibrium radius of gyration and therefore the above two alternatives are equivalent. In 
our previous study of linear chains under shear flow we observed apparent scaling exponents 
I'a for the eigenvalues of the gyration tensor with (A^ — 1) which were different from the 
equilibrium one in the flow plane directions while in the neutral (out of plane) direction the 



equilibrium scaling was preserved within error bars [^^jbH- present study on longer 

chains we find that at (3 = 3.2 the equilibrium scaling Ga oc A^^^ is preserved in the neutral 
direction [a = III) and in the in-plane compression direction [a = II) (see figure In 
the elongation direction an apparent scaling exponent higher than the equilibrium value is 
still observed [uj ~ 0.7) although the noise is considerably higher than in the other two 
principal directions. The results for the chain without EV interactions mentioned above 
demonstrate that such apparent scaling is related to finite extensibility effects rather than 
to EV interactions. For this reason we investigate universality in terms of qRg. 

In figure ^,b,c we plot, at p = 3.2 and for 30 < A^ < 300, {qRg/ \f?,)^/^S{q)/N versus 
qRg/ \/3 for a = III, II, I respectively and we compare with the MD data for N=50 at 



the same P from ref. [38|. For the latter we have used the observed equilibrium exponent 



t^MD = 0.57 rather then the Flory classical value up = 0.6. Note that in that model 
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the solvent was explicitly considered so that deviations from the mean field value are not 
contradictory. We observe that, with the exception of = 30 in the neutral direction, 
universality is indeed observed within error bars in all directions. A general feature of those 
curves is the presence of a quite slow crossover back to the equilibrium behavior at large g's 
which would be represented by an horizontal line. That = 300 is still too short to cover 
all length scales of the problem is clearly shown by the fact that even for this length the 
equilibrium behavior is not completely recovered before the chain finite extensibility shows 
up, a sign that finite chain length effects are still present. While in the neutral direction the 
effect of the flow is tiny, in the in-plane directions the flow strongly changes the scattering 
function. We note that the noise level is considerable higher in the elongation direction than 
in the other directions. Moreover statistical convergence of the results at low {3 appears 
to be considerably slower than at high (3 and at equilibrium (see also table PTI] ). This is 
compatible with the possible appearance of new characteristic relaxation times in weak flow 
as indicated by recent experiments [pOf . 

Finally we note that the structure factor obtained in our previous MD simulation on a 



different chain model embedded in an explicit atomic solvent under shear flow |3^, follows 
very closely in all directions the universal behavior predicted by the present BD simulations. 
As the former study includes automatically the hydrodynamics interactions [0,0 while 
the present calculation do neglect them from the start, it appears that HI play a secondary 
role on the chain structure under shear flow, once one works at a flxed value of the suitably 
deflned reduced shear rate. Because of the hmited length of the largest chain studied by 
molecular dynamics (N=50), we could not detect on the chain structure factors the slow 
trend back to the equilibrium behavior at large q's, which explains our previous analysis in 
terms of anisotropic scaling laws. 

In order to test the "shear" blob hypothesis we report in flgures P|P,pil|, {qRg)'^S{q)/3N 
versus qRg/ for A^ = 100 and A^ = 300 at P = 3.2 and for A^ = 100 at equilibrium, and 
we compare with the ideal chain behavior. In all directions Rouse behavior describes the 
excluded volume chain up to qRg/ y/S ~ 3. Similarly to the stretched chain case 



we 



clearly see in the compression and the neutral directions a crossover from ideal statistics 
at low g's to EV statistics at high g's. The inset of EV statistics is at qRg/ a/3 ~ 9 in the 
two directions. Note that for A^ = 100 the EV regime is completely missing and even with 
A^ = 300 it is very narrow in the compression direction. In the elongation direction the noise 
level is quite high in the Rouse (small q) regime. As stated above, in this direction flow 
effects extends to higher g's than in the other directions, above qRg/^/S = 9. On the other 
hand, the finite extensibility takes place around qRg/ ^/S = 20 for A^ = 300 and therefore 
the EV regime cannot be observed in this direction for such chain length. 

Onuki's phenomenological model is based on the assumption that below a charac- 
teristic length scale ^ the effects of the flow are negligible. The length scale separation is 
dictated by the typical relaxation time of chain density fluctuations at that scale. At a 
given shear rate 7, the crossover occurs at a scale ^ = An'/, such that the longest relaxation 
time of the blob is Tc = C,'^/{6Dc) = where = Dq/ti/. is the diffusion of the blob 
as a whole and v' is the dynamic scaling exponent. Working at flxed reduced shear rate 
(3 = TNi-i the number of monomers per blob results in Uc = N/ (3^^^"^"'^'^ ^ so that the number 
of blobs per chain is = N/n^ = /?V(2v+v') ^ ^ ^pu/{2u+u') ^ RJ f]^/i-i^+^') , For 

our model {A = 0.40448, = 3/5, u' = 1) at /3 = 3.2 we have = 1.7, Uc = 177 and 
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^ = 9.02 for = 300. The crossover from ideal to EV statistics in figure ^ is observed 
around qRgl\f?) = 4.0 which for = 300 provides Qc = 0.56 ~ 5/^ ~ Svr/^. Therefore 
the shear blob hypothesis is compatible with our results although a more firm test on the 
basis of longer chains studied at various values of f3 is still missing. To finally prove that 
our chains are too short for this purpose we show in figure |Tl| the Kratky plots in the three 
principal directions for = 300 and /3 = 3.2 and (3 = 10.0. It is evident that ai j3 = 10.0 the 
finite extensibility effect interferes with the EV regime giving rise to anomalous behaviours. 



VI. CONCLUSIONS 

In order to study the structure of a single chain subjected to a homogeneous and steady 
shear flow and to elucidate the effects of excluded volume interactions we have computed 
the scattering function of the continuous Rouse chain model in flow and we have performed 
BD simulations of a chain model of Fraenkel springs with fully developed EV interactions. 

Concerning global properties, our simulations confirm that at a given reduced shear rate, 
EV chains are less deformed and less oriented along the fiow lines than chains at the B point 

We have shown that, even in the ideal case, the chain structure in shear fiow cannot 
be explained by the known behavior of a chain stretched at its ends. In particular we 
have shown that the structure factor of the Rouse model under SF for the scattering vector 
lying along the in plane elongation direction is very different from the one of the Rouse chain 
stretched at its ends with q along the stretching direction. It does not present the oscillations 
typical of permanent extension, but it is a monotonous curve qualitatively similar to the 



equilibrium Debye curve (see figure p!0| ). In the elongation direction, the onset of the well 
known behavior signaling the static scaling of ideal chains is shifted to higher g's than 
at equilibrium and is strongly dependent on the shear rate. In the compression direction it 
is at the same location as at equilibrium and it is almost insensible to shear rate. 

Adding the EV interactions, our results suggest that the universality of ^(q, A^, 7) in 
terms of x = q-R^ and /5 only {Rg is the equilibrium radius of gyration and /5 is computed on 
the basis of the equilibrium characteristic relaxation time), exhibited by the Rouse model at 
any fiow intensity, is preserved by the EV interactions provided the chains are long enough. 
Studying chains between 30 and 300 monomers we indeed observe that the eigenvalues of the 
gyration tensor in the compression (II) and in the neutral (III) directions follow scaling laws 
with N with the equilibrium Flory exponent. In the elongation direction a larger apparent 
scaling exponent is instead measured (around 0.7), in agreement with previous MD data 
PS|^D[ which led us to propose an anisotropic scaling picture for intermediate reduced shear 



rates [^]. From the results for the same chain model without EV interactions we infer that 
this anomalous scaling arises from finite extensibility effects rather than EV interactions so 
as the equilibrium scaling should be recovered even in the elongation direction for longer 
chains. Unfortunately the dynamical character of the BD algorithm and the unfavorable 
scaling with N of the chain relaxation time does not allow us to study chains longer than 
N = 300. 

The general agreement found between our present data (where HI are absent from the 
model) and the previous non-equilibrium MD data for short chains (where HI are automat- 
ically included) strongly suggests that HI effects play a minor role on the chain structure 
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under shear flow, provided we compare results at the same reduced shear rate, which incorpo- 
rates imphcitly the HI effect through the longest relaxation time of the chain at equilibrium. 
The EV and the finite extensibility effects are by far the most relevant effects to be taken 
into account. 

For our longest chain a.t (3 = 3.2, the structure factor in the neutral (III) and in the 
compression (II) directions exhibit a clear crossover from the ideal behavior at small g's, 
well predicted by the Rouse model, to the isotropic equilibrium EV behavior at high g's, 
characterized by a power law S{q) ~ q~^^'^ with u = i^eq = 3/5. This behavior is usually 
related to the presence of the so called "blobs" , in the present case "shear blobs" . Blobs were 
detected in the chain structure factor either experimentally for semi-dilute solutions M] or 



by Monte Carlo simulations for a single EV chain stretched at its ends . The presence of 



blobs in a chain under shear fiow has been predicted theoretically but had never been 
observed so far, neither by SANS nor by simulation. The results of our BD simulations 
are compatible with the existence of blobs (see section |V B| ) but, because of finite chain size 
effects, a quantitative test of this hypothesis would require chain lengths larger than N=300. 
It is interesting to note that the largest chain considered here has a size which corresponds 
to a polystyrene chain of the order of 300 Kuhn segments of about ISA, which yields a 
molecular mass of 300.000 Daltons PB[ |. 



The present analysis should be tested by scattering experiments on dilute solutions for 
chains in good solvent subjected to steady shear flow. It would be interesting to study 
at low qRg, either by Light or by Neutron scattering, whether the Rouse model indeed 
describes the global orientation and deformation of chains for large deformations, i.e. at 
j3 » 1. Remaining discrepancies could be interpreted as a direct measure of FE and 
HI influence on these structural properties far from the equilibrium structure. Ultimately, 
the most spectacular experimental result we seek for remains a SANS study of a suitable 
"polymer-good solvent" pair which would allow us to follow the form factor of a polymer of 
a few million Daltons under shear flow at large (3 over a broad q regime. This would allow a 
direct observation of the crossover from ideal chain to EV chain statistics across a crossover 
diffusion vector q^ characterizing the shear blob dimension in a particular direction. 

Concerning future simulation work, a more quantitative test of HI effects on the chain 
structure is in order and could be performed by BD on short chains. 
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APPENDIX A: IDEAL CHAIN UNDER STEADY SHEAR FLOW 

In this appendix we derive the steady state distribution function of the continuous gaus- 
sian chain model under an homogenous shear flow, eqs. (^, (0) and (||). 

The dynamics of the model under a generic flow fleld, specifled by the velocity gradient 
tensor k, is deflned by the Langevin equation 
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dt 

with the boundary conditions 



e^ = X^^ + f.(t)+ek(t).R„(t) 



dn? 



'dRnity 




'dRnity 




dn 


n=0 


dn 


n=N 



f„(t) is the gaussian random force acting on the chain and is defined by 

< Ut) > = 

< f„(t)f„(t') > = 2kBT5{t - t')5{n - m)l 



(Al) 



(A2) 



(A3) 
(A4) 



where < ... > indicates an average over the distribution of the random noise and 1 is the 
unit tensor. As it is well known the above Langevin dynamics corresponds to the following 
diffusion equation (Smoluchowsky) for the distribution function in the chain configurational 
space [|,^ 



= - dnV 



N 



(9^ 
'di 

with the dynamics of monomers given by 



knT 



^ din? ^ 
Let us introduce now the normal modes defined as 



Xp(t) 



N Jo 



N 



dn cos [-j^j ttnit) 



p = 0,l,2, 



in terms of which the chain positions are given as 
R„(t) = Xo(t) + 2 f; CO. (^) X,{t) 

The completeness relations for the normal modes basis are 



0<n< N 



(A5) 



(A6) 



(A7) 



(A8) 



Jpq 



N Jo 



N 



dn cos 



N6{n -m) = 1 + 2 J2 cos 



/ p7m\ 



f qnn 



\ N 



cos 



N 
/ P'Km\ 



) 



(A9) 
(AlO) 



Note that, for homogeneous fiows which we consider here, the normal modes of the dynamics 
are the same as at equilibrium. 

The Langevin equations for the normal modes are 



^pdtXpit) = -XpXp(t) + fp(t) + ^pk ■ Xp(t) 
with the following definition of the symbols 



(All) 
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Xp 



p= 1,2,3... 



< fp(t) > = 0, < fp(t)f,(t') >= 2^pkBT6pg6{t - t')l 
The distribution function of normal modes is factorized 



^({x},t) = n^p(Xp,t) 

p=0 

and one can easily derive the following diffusion equation for the p — th mode {p > 0) 



(A12) 

(A13) 
(AM) 

(A15) 



dt^Pp{Xp,t) = -Vx, 



(k • X, - fx,)ij, - ^Vx,i 



(A16) 



The center of mass motion is instead described by the diffusion equation for the zero mode 



k ■ Xo^o 



2e 



-VxoV'O 



(A17) 



The solution of eqs. (|A16|) is a multivariate guassian 



1 



1/2 



e 2 



^(27r)3rfet[ap(t)]^ 

where the time dependent covariance matrix ap{t) is given by (see eq.(7.161) in ref. 



(A18) 



ap(t) 



2p- 



dt'^^E{t,t') ■ E^(t,t')e"^ 



(t-f) 



(A19) 



Here ti = ^i/xi is the relaxation time of the first normal mode and the tensor E(t, t') is the 
deformation tensor of the external field defined by 

E(t,t') = 1+ t dt'\{t") (A20) 
Jt' 

For the steady shear flow to which we are interested in this work we have 



/O 7 0\ 


Vo 0/ 



/I i{t-t') 

E(t,t')= 1 
\0 1 

so that the covariance matrices Up results time independent and takes the form 



(A21) 



(A22) 
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Xp 



k - 



where we have defined the matrix f3p through the second equivalence (see eq.( 
that 



det{ap 



7 



2p2 



(A23) 
Noting 

(A24) 



the distribution function of ihep — th normal mode under steady shear fiow results in eq.(0). 
It is easy to check that eq.(0) reduces to the standard equilibrium distribution for 7 = 0. 



APPENDIX B: DERIVATION OF THE SUMS 



In this appendix we outline the method to compute the sums in eqs. (|l^) 



Sk{a,h) = ^^^[cos{p-Ka) ~ cos{pTih)f fc = 0,1,2 (Bl) 



p=i 



where < (a, 6) < 1. It is evident that Sk{a,b) = Sk{b,a) > and the equality holds for 
a = b only. Developing the square inside the sum in eq. (pi]) and using standard relations 
among trigonometric functions we obtain 

Sk{a, b) = (Tfc(O) + ^ [(Tfc(a) + ak{b)] - a,. {^-^ - ^fc (^2) 

with 

= (B3) 
p=i y 

so that the problem reduces to compute (Tfc(x) for A; = 0, 1, 2. Note that ak{x) is an even and 
periodic function of x with unitary period. Eq. ([B3D expresses the Fourier series a function 



fk{x) with these properties, after removing the p = term. Therefore we have to seek for a 
function fk{x) that satisfies the following conditions 

2p2+2fe r dxfk{x) cos{2npx) = 1 p>l (B4) 

Jo 

1 

dxfk{x) sin{27ipx) = 







Applying repeatedly integration by parts we easily obtain the following identities 

/ dxfkix) cos{2TTpx) = — — — / dxfKx) d^cos{2TTpx) (B5) 
JO 27rn r jo 



(27rp)2 Jo 
1 

(27rp)4 
1 



(27rp)6 Jo 



f dxf'^'ix) d^cos{2'Kpx) (B6) 

JO 

dxfl{x) dxCos{2T[px) 
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having imposed the boundary conditions (/fc(0) = for the function and for its deriva- 

tives up to order five. 

Let us now consider the three k values separately. For k = we have 



1 




1 



27r2 Jo 
1 



(27rp)2 Jo 
/o(0) = /o(l) 



(ixf^i^x) dxCos{2npx) 
dxfQ{x) dxSin{2TTpx) 



(B7) 



p>l 



where we have used the first of the equalities in eq. (|B5|) and a similar one for the sinus 
function. These relations can be satisfied by choosing fo{x) = Cq + Ci x + C2 with 
Co = 0, Ci = — TT^ and C2 = vr^ which provides /o(x) = 7t'^{x'^ — x). The function cro(x) is then 



aQ[x) 



1 1 
dx /o(x) = 7r^(x^ — X — -] 



and, after eq. (p2| 



So{a,b) 

which is the wanted result for k = 0. 
For k = 1 we must solve the problem 

1 '•I 



TT 



\a-b\ 



(B8) 



(B9) 



1 




in^ Jo 



(27rp)4 Jo 
/i(0) = /i(l) 



dxfi{x) dxCos{2'Kpx) 
dxfi{x) dxSin(2TTpx) 



(BIO) 



p > 1 



where the second equalities in eq. (|B5|) has been used. These relations can be satisfied 
choosing /f (x) =03 + 04^ with C4 = — Svr^ and C3 = 47r^. Integrating three times over x and 



imposing the boundary conditions on /( and /i we obtain fi{x) = 7r'*(2 x^ 
which provides 

= - / dx fi{x) 
Jo 

Using eq . ( p^) we finally obtain 



X 



TT 



4 1 



3 '30 



x^ + 2x^ - x^) 



Si{a, h) = —{a — b) 



(a + fe) 1 



a + 



\a — b\ 



In the case k = 2 the relation to be satisfied are 

2 



1 







(27r)6 Jo 
1 



(27rp)6 Jo 
/2(0) = /2(1) 



dxfl{x) dxCos{2'Kpx) 
1 

dxfl{x) dxSin{2TTpx) 



x^)/3 



(Bll) 



(B12) 



(B13) 



p > 1 



19 



A procedure similar to the previous cases provides f2{x) = 7r^x'^(2x'^ - 6x^ + 5x^ - l)/45 
and therefore 



cr2{x) 



71 

45 



5 



+ x^{2x^ -6x^ + 5x^ - 1) 



(B14) 



and 



S2{a,b) = ^{a-bf{\a-bf + 5{a + b) \{a + b){2 + a'^ + b^) - Sa^ - 2ab - 3b^]} (B15) 
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TABLES 



TABLE I. Equilibrium simulation details and results, tn is defined through the scaling laws 
given in the text. Nt is the total length of the simulation run in time steps. T/tn is the total 
simulation time over the longest chain relaxation time. 



N 


AtxIO"^ 


Rg 


R 


TN 


T/tn 


9 


2 


1.413(5) 


3.44(1) 


3.04 


165 


20 


10 


2.340(2) 


5.83(1) 


19.1 


131 


30 


50 


3.07(2) 


7.77(5) 


47.6 


263 


50 


50 


4.16(1) 


10.4(1) 


149 


84 


100 


300 


6.410(1) 


16.27(2) 


677 


111 
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TABLE II. Flow simulations details and results for the stiff spring chain model without EV 
interactions. N = Rouse means data for the continuous model of section |l|. T/tn is the total 
time of the simulation over the equilibrium longest relaxation time of the chain 



N 






mo 




6Gi 




SGni 


100 


3.2 


200 


2.8(3) 


1.7(1) 


2.6(3) 


-0.33(1) 


- 0.09(1) 


300 


3.2 


120 


2.7(3) 


1.78(4) 


5.4(1.2) 


-0.32(2) 


- 0.02(1) 


Rouse 


3.2 




2.5 


1.75 


5.00828 


-0.32714 





100 


10.0 


200 


3.7(1) 


2.25(3) 


21(2) 


-0.44(1) 


-0.18(1) 


300 


10.0 


60 


2.8(2) 


1.97(6) 


35(3) 


-0.41(2) 


-0.03(2) 


Rouse 


10.0 




2.5 


1.75 


46.06164 


-0.34736 






TABLE III. Flow simulations details and results for the stiff spring chain model with EV 
interactions. Symbols are like in table II. Data from ref.[38] are also reported for comparison. 



N 


/? 


T/tn 


mo 


■mo 


6Gi 


SGn 


6G111 


20 


3.2 


131 


4.2(1) 


2.60(6) 


1.41(3) 


-0.430(6) 


-0.0061(2) 


30 


3.2 


263 


4.0(1) 


2.5(2) 


1.5(1) 


-0.42(1) 


-0.16(2) 


50 


3.2 


168 


3.7(4) 


2.5(2) 


1.8(2) 


-0.40(3) 


-0.17(1) 


100 


3.2 


111 


4.1(6) 


2.3(3) 


1.4(1) 


-0.40(1) 


-0.21(2) 


300 


3.2 


98 


3.3(3) 


2.1(3) 


2.5(3) 


-0.47(2) 


-0.14(2) 


50^^ 


3.2 


73 


4.6(7) 


2.1(9) 


1.8(2) 


-0.47(1) 


-0.12(1) 


100 


10.0 


75 


6.2(3) 


3.9(1) 


5.2(2) 


-0.60(4) 


-0.34(5) 


300 


10.0 


59 


5.2(2) 


3.4(1) 


8.1(8) 


-0.57(1) 


-0.24(2) 


30^ 


10.0 


167 


6.2(4) 


3.4(3) 


5.3(3) 


-0.62(3) 


-0.40(4) 



'^data from ref. 
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FIGURES 




FIG. 1. Structure factor of the continuum Rouse chain model under shear at various values of 
f3 and for the scattering vectors along the principal axes of the gyration tensor in the flow plane. 
P = 0.0, thick line, /3 = 1.0 dot-dashed lines, P = 3.2 thin line, (3 = 10.0 dashed, and P = 100.0 
long dashed. Curves below the equilibrium one are in the elongation direction I, while curves above 
the equilibrium one are in the compression direction II. 
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FIG. 2. Equilibrium static scaling for the end-to-end distance (circles) and the gyration radius 
(squares). 
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FIG. 3. a) Static structure factor at equilibrium for various chain lengths from = 9 to 
N = 100. b) Kratky plot S{x)x'^ /N in terms of x = qRg. The tick dashed curve is the ideal 
(Debye) behaviour. The dot-dashed straight line represents the x^^'^ behavior expected for EV 
chains. 
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FIG. 4. Scaling of the internal dynamics provided by the dynamics structure factor as explained 
in the text. A/" = 20 circles; N = 30 squares; N = 100 diamonds. Closed symbols are for I = 0.1 
while open symbols corresponds to / = 0.6. 
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FIG. 5. Orientation angle (closed symbols) and extinction angle (open symbols) at /3 
versus N. 
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FIG. 6. Scaling of the eigenvalues of the gyration tensor in the molecular reference frame. The 
equilibrium behavior is represented as a thick dashed line. Continuous lines arc power law fits to 
the data. The apparent scaling exponents in the three directions are reported on the plots 
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FIG. 7. Structure factor for the EV model at /? = 3.2 and for various N as indicated in the 
figure. MD data from ref.[38] for iV = 50 at /? = 3.2 are reported as a thick dashed curve. 
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FIG. 8. Ideal Kratky plot in the neutral (out of plane) direction III at /3 = 3.2 for N = 100 and 
N = 300. The thick dashed curve is the equilibrium behavior for N = 100. The Rouse behavior 
and the asymptotic EV scaling are represented as continuous lines. 
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